Wavefunctions and counting formulas for quasiholes 
of clustered quantum Hall states on a sphere 



O 
O 

(N 

>>: 



i 

C 

o 
o 



> 
00 

o 
o 

-I— > 



i 

G 

O 

o 



N. Read 

Department of Physics, Yale University, P.O. Box 208120, New Haven, CT 06520-8120 

(May 5, 2006) 

The quasiholes of the Read-Rezayi clustered quantum Hall states are considered, for any number 
of particles and quasiholes on a sphere, and for any degree k of clustering. A set of trial wavefunc- 
tions, that are zero-energy eigenstates of a k + 1-body interaction, and so are symmetric polynomials 
that vanish when any k + 1 particle coordinates are equal, is obtained explicitly and proved to be 
both complete and linearly independent. Formulas for the number of states are obtained, without 
the use of (but in agreement with) conformal field theory, and extended to give the number of states 
for each angular momentum. An interesting recursive structure emerges in the states that relates 
those for k to those for k — 1. It is pointed out that the same numbers of zero-energy states can 
be proved to occur in certain one- dimensional models that have recently been obtained as limits 
of the two-dimensional k + 1-body interaction Hamiltonians, using results from the combinatorial 
literature. 



Non-Abelian quantum Hall states have been a subject 
of great interest for some time [1] . Moore and the author- 
showed that non-Abelian statistics, in which adiabatic 
exchange of well-separated quasiparticles produces a ma- 
trix operation on a space of degenerate states, was a pos- 
sibility in a condensed matter system, and gave an exam- 
ple, now often termed the Moore-Read (MR) state. They 
related the deep structure underlying the statistics and 
other aspects of the quasiparticles of these phases to con- 
formal field theory, and to what are now called modular 
tensor categories. Many trial wavefunctions of quantum 
Hall states can be related to conformal blocks (i.e. chiral 
parts of correlation functions) in the corresponding con- 
formal field theory. A further sequence of examples was 
constructed by Rezayi and the author (RR) [2]. These 
are parametrized by two integers, k > 1 and M > 0; the 
filling factor for each is v — k/(Mk+2), and they describe 
bosons for M even, fermions for M odd (the case k = 1 
is the Laughlin state, while k — 2 is the MR state). For 
M = or 1, trial ground state and quasihole states can 
be obtained as zero-energy eigenstates of a Hamiltonian 
that consists of a purely k + 1-body interaction, acting 
within the lowest Landau level. (For M > 0, the trial 
wavefunctions are those for M = multiplied by the fac- 
tor Yli<j( z i~ z j) M > wner e Zj = Xi+iyi is the complex co- 
ordinate representing position (xi,yi) of particle i in two 
dimensions. These trial wavefunctions are zero-energy 
eigenstates for a Hamiltonian that contains two-body as 
well as k + 1 body interactions [2].) For M = 0, the 
wavefunctions are closely related to conformal blocks of 
SU(2) level k current algebra conformal field theories [2]. 
While the explicit trial ground state wavefunctions were 
obtained, for k > 3 the trial wavefunctions for quasiholes 
were not, except for some of those for small numbers of 
quasiholes [2,3]. The number of degenerate zero-energy 
states for fixed positions of the quasiholes is basic to the 
theory of non-Abelian statistics, and for the MR and RR 
states these degeneracies have been studied numerically 



and using conformal field theory methods [4,5,2,6-8]. 

The MR state is now believed theoretically to describe 
the phase observed at filling factor v = 5/2 in electron 
systems [9], while there is evidence that the RR k = 3, 
M = 1 state occurs at v = 12/5 [2]. For bosons, there 
is numerical evidence that the sequence of RR states is 
relevant to rotating cold, trapped, bosonic atoms [10,11]. 
Interest in non-Abelian statistics is presently being fueled 
by the possibility of topological quantum computation 
[12]. 

In this paper, we return to the question of finding and 
counting the zero-energy states of the k + 1-body inter- 
action for bosons, for general k. We work in the lowest 
Landau level on the sphere, with quanta of magnetic 
flux piercing it [13]. The approach is a generalization of 
one used to prove linear independence of a proposed set 
of functions for the k = 2 case in Ref. [5] , and is very close 
to work [14] on the same counting problem for the limit 
N$ — > oo, which has been reviewed, extended, and ap- 
plied in Ref. [15]. (The authors of the first paper in Ref. 
[14] discovered the relation of the symmetric functions 
that vanish when k + 1 coordinates coincide to conformal 
blocks of SU(2) current algebra at level k, which was dis- 
cussed independently in Ref. [2] . From the quantum Hall 
point of view, the results on the — > oo case describe 
edge excitations.) Our main results are: (i) a complete 
and linearly-independent set of explicit wavefunctions for 
zero-energy states of the k + 1-body interaction for any 
number of particles N and any number of flux quanta N^, 
or alternatively, any number n of quasiholes, as defined 
by 



2 n 
N$ = -N + M(N-1) + - 



(1) 



(ii) explicit generating functions that encode the number 
of zero-energy states for any N, n, decomposed accord- 
ing to orbital angular momentum L z (and hence L also); 

(iii) an interesting recursive structure emerges, that re- 
lates the quasihole problem for k to the similar problem 
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for fc— 1 at a reduced number of particles and flux quanta, 
and might be useful in making further progress in under- 
standing the quasihole states. The arguments are short 
and self-contained, and make no use of conformal field 
theory. In an appendix, we also consider quasihole wave- 
functions of a natural form that generalizes those in Ref. 
[3], but for these our results are less conclusive. 

We want to find the complete set of zero-energy states 
for the k + 1-body Hamiltonians in Ref. [2] . The cases 
with M > can be found by multiplying those for M = 
by the Laughlin-Jastrow factor Yii<j ( z i~ z j) to the power 
M, so we can focus on the M = case (hence, on bosons). 
Then the problem reduces to the following: find all the 
symmetric polynomials in N complex variables Zi that 
vanish when any k+ 1 of the Zi are equal. For the present, 
working on the sphere means that the polynomial is of 
degree at most in each Zj [13]. We will denote the 
vector space of symmetric polynomials with these prop- 
erties by VN,N#,k an d a typical member of this space by 
^N,N^,,k (the polynomials 4> can be related to wavefunc- 
tions \I> on the sphere or plane as explained in Ref. [5]). 

Consider a function ^N,Nj„k- While it vanishes when- 
ever k + 1 particles coincide, it may not vanish when 
only k coincide. Let z\, z 2 , Zk tend to the same 
value, which we will denote by Z\. If it vanishes in this 
limit, then (since it is totally symmetric) it vanishes when 
any k coordinates coincide, and belongs to Vjv.iv^.fc-i- If 



The notation and the justification for each of the factors 
is as follows. The first product is over pairs of cluster co- 
ordinates Z p , labeled by p, p', . . . , = 1, 2, . . . , (N — F)/k 
(the number of clusters). The function vanishes when 
any k + 1 particles coincide, so it must vanish when any 
two clusters coincide, and since it is symmetric under ex- 
changes of clusters, it must contain a factor (Z p — Z p >) 2 
for each pair of clusters. In fact, because the original 
wavefunctions are symmetric functions for particles on 
the sphere, it must vanish as the 2fcth power; the proof 
of this is given in an appendix to this paper. The last 
factor ^ F,Nj,-2(N-F)/k,k-i i s a function only of the F 
"unclustered" coordinates, so is independent of the clus- 
ter coordinates. By assumption, it vanishes if any k of the 
remaining F coordinates coincide, so can be any function 
in Vp^N^-2(N-F)/k.k~i- The residue must also vanish if 
any of the unclustered particles approaches a cluster, so 
it must contain Z p — zx-F+m for the pth cluster and the 
mth unclustered particle (m = 1, 2, . . . , F). In fact, it 
must vanish quadratically, which again is proved in the 
Appendix. This form implies that the maximum degree 
in the unclustered coordinates in the remaining factor is 
— (N — F)/k, as recorded in the second subscript on 



it does not vanish, then its value (corresponding to the 
"residue" in Ref. [5]) is a function of Z\ and of Zk+i, ■ • ■ , 
zn, and is symmetric in the latter. In this residue we 
may let another set of k coordinates, say Zk+i, ■ ■ ■ , z 2 k 
coincide at Z 2 . If the function vanishes in this second 
limit, then the residue consists of a function of Z\ times 
a member of Viq-k.N^.k-i- If it does not vanish in the 
second limit, we obtain a second residue, to which the 
procedure can be applied again. 

In this way the space of functions is divided into a 
nested sequence of subspaces within subspaces (called a 
"filtration" in Ref. [15]). Each such subspace contains 
all the wavefunctions that have a non-vanishing residue, 
which is a function of Z\ , Z 2 , Z%, . . . , (each of which 
are the locations of (N — F)/k "clusters" of k particles), 
and of the remaining F coordinates zn—f+i, ■ ■ ■ > z n of 
the F "unclustered" particles; these residues vanish if 
any k of the latter set coincide. Each subspace also con- 
tains functions with vanishing residues, which are also 
members of a subspace with fewer clusters. Thus, the 
(N - F)/kth subspace [(N - F)/k = 0, 1, . . .] can be 
defined as the space of functions that vanish identically if 
(N — F)/k + 1 clusters of k particles each coincide at Z\ , 
. . . , Z(jv-F)/k+ij f° r au values of these coordinates and of 
the remaining z's. We will show that the non-vanishing 
residues of functions in this subspace must take the fol- 
lowing form: 

(2) 

I 

^F,N^,-2(N-F)/k,k-i- Finally, the residue can contain a 
further factor of a symmetric polynomial in the cluster 
coordinates. As the degree in any of the cluster coordi- 
nates is at most kN^,, it works out that the maximum 
degree in each cluster coordinate in this symmetric poly- 
nomial factor is n. Each quasihole effectively contains 
1/k of a flux quantum or vortex [2]. In general, n need 
not be divisible by k (in Ref. [2] , only the special cases N 
and n divisible by k were considered). For to be an in- 
teger, we require that 2N + n be divisible by k (hence for 
k even, n must be even), and dimVjv i jVj,,fe (the dimension 
of Vft^jv^fc) i s nonzero only when n > 0. The symmet- 
ric polynomials in the cluster coordinates of this degree 
can be expressed using the product rj „ ;(^ p » — w{) as a 
generating function (by expanding it in symmetric poly- 
nomials in the variables wi). These are then "coherent 
states", in which the coordinates wi can be interpreted 
as the positions of the quasiholcs. The function vanishes 
when any k particles approach awi, as it should [16], and 
this can be taken as a definition of the quasihole coor- 
dinates for these quasiholes, which have charge —1/2 (in 
units of particle number). 

Let us suppose that the problem has already been 



I 

^N.N^,k —> Y\_ ( Z P ~ Z P') 2k Y\{ Z P" - w l) ] I ( Z P"' - z N-F+m) 2x &F,N<p-2(N-F)/k,k-l- 
p<p' p" \l p'" ,m 
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solved for values k' = k — 1. Then we have a set of 
linearly independent functions <J/ F,N 4 ,-2(N-F)/k,k-i that 
span V FyN<j> _ 2 (N-F)/k,k-i- Then by expanding in the 
quasihole coordinates u>/, we have a set of non- vanishing 
residues that are functions of the cluster coordinates and 
the unclustered particle coordinates, and clearly are lin- 
early independent. This implies that if for each of these 
residues, there is a wavefunction from which it arose by 
taking limits, then these wavefunctions (which are in 
VN,N^,k) are linearly independent. The existence of a 
wavefunction that reduces to each non- vanishing residue 
in the limit will be shown in a moment. 

Completeness of the set of functions can also be proved 
by induction. As the form of the non-vanishing residues 
has been fixed, we must show that there is a one-one 
mapping from the wavefunctions to the non-vanishing 
residues. This is done by induction on the number of 
clusters. For the residues with no clusters (F = N), 
the residue reduces to ^ N^^.k-i, for which a complete 
set is assumed to have already been found. Then for 
the residues with one cluster (F = N — k), if there are 
two wavefunctions with the same residue, then the differ- 
ence of these functions is a wavefunction with no clusters 
(since it vanishes when k particles coincide), and these 
are also in the space Vjv,iVj,,fc- This argument can be re- 
peated, and so by induction we find that there is (up to 
addition of a wavefunction for fewer clusters) at most one 
wavefunction for each non-zero residue, that is (assum- 
ing existence) the dimension of the space of non- vanishing 
residues and that of the space of wavefunctions are the 
same for each N, N^, and k. 

To prove the existence of a zero-energy wavefunction 



that tends to each non- vanishing residue in the limit, we 
construct them. We use the following labeling for the par- 
ticles: The particle number N is partitioned into m a > 
clusters of size a = 1, . . . , k, that is N — am a (there 
will be a set of wavefunctions for each such partition). 
One particle is assigned to each box in a Ferrers- Young 
diagram in which there are at most k columns, and which 
consists of m a rows of length a for each a, arranged as 
usual with the row lengths weakly decreasing as one goes 
down the diagram. Then (following Rcf. [15]) the parti- 
cle coordinates will be written where i, j label rows 
and columns, respectively (as for a matrix, i increases 
down the columns, and j increases to the right along the 
rows) of the rectangular block consisting of all rows of 
length a = 1, . . . , k. Thus for each a, j = 1, . . . , a, 
and i = 1 , . . . , m a . The rows of the diagram are ordered 
with (a, i) above (/3, i'), written (a,i) > ((3,i'), if either 
a > (3, or a = [3 and i < i' (note the direction of the last 
inequality) . The wavefunctions will be written down in a 
form that deals with all levels of the recursion at once. It 
will turn out that (N — F)/k = , while the remaining 
m a represent, similarly, the numbers of clusters of sizes 
a smaller than k. Similarly, the wavefunctions involve 
not only quasihole coordinates (I = 1, . . . , = n), 
which can be identified with the wi in the residues, but 
also further coordinates w\°^ for all a = 1, . . . , k — 1, 
which play a corresponding role in the parametrization 
of the wavefunctions for smaller k that appear in the re- 
cursion. Here for each a the range is I = 1, . . . , n a , with 
n\ < n 2 < ■ ■ ■ < n k . 

The wavefunctions are 



I 

* ...,z N -, {«,«>, , =5,( n n (*s ) - *$)(*sii - nnn (aj) (^ - ^) 

{(aA)>(f3A>)j=l,...,f3 a i,j I 

(3) 

In this expression, S z is the symmetrizer over all the TV particles, and z\°^ +1 — z^f . In each product nl"'^; ' runs 
over 

I = rij-t + (1 - (5^1 ) ^2 mp + l,...,rij+ mp (4) 

P:j-l<P<a P:j-l<P<a 



(where n = 0). To ensure that each function is of the 
same degree = (2N + n)/k — 2 in all coordinates Zi, 
we require 

2 ^2 mp + n a - n Q _i = (2N + n)/k (5) 

p:p>a 

for all a = 1, . . . , k. Also, for each a, the wavefunctions 
are of the same degree m a in w\°^ for all I = 1, . . . , n a . 

The structure of the z — z' factors inside the sym- 
metrizer is based on that in Ref. [15], and for the factors 



that connect two rows (a,i), ((3,i') with a — (3 is the 
same as in the ground states (i.e. the case = N/k, 
m a = for a < k) in Rcf. [2]. These particle-particle 
factors are of two types: one type connects two particles 
(and vanishes when their coordinates coincide) whenever 
they are in the same column in the diagram. The other 
type connects two particles in different rows only if the 
one in the higher row is one column to the right of the 
lower one; when the rows are of the same length, this is 
interpreted cyclically along the rows. The following ar- 
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guments are a version of the arguments for existence of 
a wavefunction mapping to each residue in Ref. [15], and 
generalize an argument for the ground state wavefunc- 
tions in Ref. [2]. 

The wavefunctions in eq. (3) vanish when any k + 1 
coordinates are equal, because in each term inside the 
symmetrizer, at least two must be in the same column of 
the diagram (i.e. have the same value of j). They also 
vanish if k particles coincide at Z\ say, unless the par- 
ticles lie in the same row (which must be of length k). 
First, these particles must be in distinct columns for the 
wavefunction to have a possibility of not vanishing, and 
hence one must be in column k, and so in the first rrik 
rows. Then in order for the function not to vanish as the 
particles come together, the one in column k — 1 must be 
in the same or a higher row, and so on, until in the first 
column the particle will be in the same or a higher row as 
that in column k. If the row is higher, the wavefunction 
will vanish because of the cyclic connection of the first 
and last columns. Hence all k must lie in the same row. 
Repeating this argument as another k particles tend to 
Z2, and so on, one finds that in the clustered limit, when 
there are Wfc clusters with k coinciding coordinates in 
each, the particles in each cluster came from the same 
row in each term that contributes a non-vanishing term 
to the residue. The residue vanishes if any k of the re- 
maining coordinates coincide, and vanishes as (Zi — Zj) 2k 
if any two clusters of k coincide. The wavefunctions also 

(k) 

vanish if k particles coincide with a w t . Again, the k 
particles must be in the same row, and then exactly one 
of them appears in a factor (z — w\ k ^), so the wavefunc- 
tion vanishes, and the residue does likewise. Finally, the 
residue vanishes quadratically as (Zi — z') 2 for any of the 
remaining unclustered coordinates z', as each of them is 
connected to exactly two particles in each row of length k. 
This establishes that the residue is of the form in eq. (2), 

with $f F N > k' being a function in the remaining z's and 

' 4, ' 

w's of the same form as in eq. (3), with parameters as in- 
dicated. Then arguing by induction on k establishes that 
there is a wavefunction for each non- vanishing residue, as 
claimed. 

Taking all the arguments together, we can conclude 
that the wavefunctions in eq. (3), though written as over- 
complete sets of "coherent states" using auxiliary coordi- 
nates w\ a>> , span a complete set of linearly- independent 
zero-energy wavefunctions. Because the residues are in- 
variant under permutations of the w^'s among them- 
selves (for fixed a), no states are lost if the wavefunc- 
tions are symmetrized over the w^'s also, by applying 
Y\ a S w ( a ) . Then the functions can be expanded as linear 
combinations of polynomials in the z's times products 
over a of symmetric polynomials in the w^'s for each 
a, and the polynomials in z for each distinct product of 
symmetric polynomials in the w's form a complete and 
linearly independent set of wavefunctions. 



For k = 1, the wavefunctions are the Laughlin quasi- 
hole wavefunctions. For k = 2 they are the same com- 
plete set of linearly-independent functions found in Ref. 
[5] (and for N even and n — 2 quasiholes, in Ref. [1]). 
In the latter, the unclustered (unpaired) particles were 
viewed as fermions occupying zero modes, but we now 
see that they can also be viewed as Laughlin quasiholc 
wavefunctions for v = 1/2 involving F particles and ri\ 
quasiholes. For general k, they agree with Ref. [2] for 
n = k, N divisible by k. 

We can now find the "number of states" #jv,Ar^,,fe = 
dim Vjv,Af^,fe, by counting the non-vanishing residues (or 
by counting the wavefunctions directly). In these first 
results, we use only the recursive structure. We expand 
the symmetric polynomials in the residues in eq. (2) in 
the «Vs, so as to obtain linearly-independent residues; 
the number of linearly-independent such functions is a 
binomial coefficient [5], and this then yields 

#N,N^,k= ( * n +U )#F^,k>, (6) 

F:F=N(modk) ^ ' 

where = iV - 2(N - F)/k = (2F + n)/k- 2, and 
k' = k — l. This has the general structure found for k = 2 
in Ref. [5]. For general fc, the structure was conjectured 
in Ref. [2], and results were obtained in [6,7] (although 
the form of the "spatial degeneracy" factors, i.e. the bi- 
nomial coefficients in this formula, seems never to have 
been derived generally). In the literature, this form is 
usually written with the notation [6,7] 

[l] = (7) 

when n + 2F is divisible by k. Remarkably, we have 
found that these coefficients, which represent the "inter- 
nal" degeneracy of the quasiholes, are themselves related 
to a problem of counting the full degeneracy of zero en- 
ergy quasihole states, but at a different value of k. This 
now allows an inductive solution for these coefficients, 
recursively in k. 

First, we can write a recursion relation for the coeffi- 
cients: 

1 J * F':F'^F(modfc<) V 11 ' 1 J *=' 

where again k! = k — 1 , and 

ri = {k'n-2F)/k. (9) 

Next we mention that for k = 1, the residues above 
are just the wavefunctions, and the function from the 
problem for k! = is simply unity. Those wavefunctions 
are the Laughlin quasihole functions. This implies that 

~ ^f,o- In this case, the statistics is Abelian, 
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and the only degeneracy is due to the different possible 
locations of the quasiholes. The recursion relation then 
yields for k — 2, 



n/2 
F 



(10) 



as found in Ref. [5]. When these numbers, or their sum 
over F, are larger than 1, they describe the degeneracy of 
states for quasiholes even when their positions are fixed, 
a characteristic of a non-Abelian state. 

When these values are substituted into the recursion 
relation for k — 3, the result can be immediately rewrit- 
ten 



{?}.- 



a,b:a+2b=F 



i-F 
3 

a 



2w-(2a+b) 
3 

b 



(11) 



which was found previously by Ardonne [7], by a much 
longer method using conformal field theory and results 
from Ref. [17]. 

The solution for arbitrary k can be found recursively, 
and parametrized in the same way as the wavefunctions, 
eq. (3), using the partitions of N and n. First, the total 
number of states can be written as 



mi ,. . .,mfc am 



n 



(12) 



Here the recursive structure is clear; notice that m k = 
(N - F)/k, n k -i = n', m fc _i = (F - F')/k', etc., so 
this agrees with preceding results. In practice we need 
the values of n a for given N and n; these are obtained 
by solving the conditions on the number of w\°^ at each 
level a in terms of the mp's: 



2N- 



a 



k 



2 i 3m ^ 

fl<a 



(13) 



for the coefficients 



, as a sum over partitions of 



for a = k, this reduces to n k = n. There is a similar form 
n 
F 

F = N — krrik, which agrees with the result of Ref. [7]. 

One further refinement is to count zero-energy eigen- 
states for each angular momentum value L z on the 
sphere. Instead of dimV^jv. k, we can choose to calcu- 
late the trace of q NN <t>/ 2 - L " ) where q is an indeterminate, 
and L z is the z-component of total angular momentum L. 
We recall that for the single-particle wavefunctions, the 
function z m contributes N<p/2 — m to L z , for m = 0, . . . , 
[13,5]. Thus (for L z eigenstates) we have arranged 
that the exponent of q is simply the total degree of the 
wavefunction in the z's. From the number of states of 
each L z , one can extract the number of quasiholc states 
of each value of L 2 also. This information can be valuable 
in numerical studies. 



Once again it is sufficient to count residues, this time of 
each degree. The total degree in the cluster coordinates 
Z p , p — 1, . . . , mfc, is km k (m k — 1) plus the total degree 
of the symmetric polynomial in the Z's produced by ex- 
panding the factor FJ ; (Z p — toj), plus the contributions 
of the factors that connect with the unclustcred parti- 
cles. The number of symmetric polynomials of degree d 
in rrik variables and of degree at most in each variable 
is given by the coefficient of the q d term in the (Gauss) 
(/-binomial coefficient (a polynomial in q), defined by 



m k + n k \ _ (m k + n k ) q \ 
n k ) {m k ) q \{n k ) q V 



(14) 



where the g-factorial (n) q \ = (n) q (n — l) q ■ ■ ■ (l) q , and 
the g-deformed integers are 

(n) q = l + q + q 2 + --- + q n - 1 = T ^-. (15) 

This combinatorial result can be understood from its own 
generating function, which is 



1 



(1 - x)(l - xq) ■ ■ ■ (1 - xq m * ) ^ \ n k 



(16) 



this identity is Heine's q-binomial formula [18]. It is eas- 
ily seen that the left-hand side of eq. (16) is the generat- 
ing function for symmetric polynomials in m k variables, 
of degree at most n k in each variable: The symmetric 
polynomials in m k variables Z p are (freely) generated 
by the elementary symmetric polynomials e\ = ^2 p Z p , 
e 2 = J2 P<P ' Z p Z p 'i ■ ■ • > e m k = Hp Z p , which are of total 
degrees 1, . . . , m kl respectively, and of degree at most 1 in 
each variable. Each factor (1 — xq m )~ l gives a series that 
contains one term x r q mr with coefficient 1, correspond- 
ing to each distinct possible factor e r m . Similar formulas 
apply to each level in the recursion, with a in place of 
k. In all these expressions, the limit q — > 1 reproduces 
earlier or well-known results, with ordinary binomials, 
factorials, etc., in place of their g-analogues. 

Since the g-binomials count the effects of the factors 
containing wj^'s, we now only require the total degree 
(summed over a) of the minimal degree polynomial in 
the Z's of all levels, that is the term of highest possible 
degree in all w's. For the given partition described by 
m a , this minimal total degree is 



am a (m a — 1) + 2am a mp, 

a Q</3 

= Y m a mf3M al 3 - N, 



(17) 



a,/3 



where M a p = min(a,/3) [14,15]. 

The result is that the trace over Vjv.jVi.fc is 



5 



tr q 



= Q 



-N 



E 



n 



m, 



(18) 



This expression counts the zero-energy states on the finite-size sphere for each L z ; furthermore, trg~ L * for M > 
is independent of M. As a test, this formula has been compared with some rows for k > 2 in the tables of numbers 
of states for each L that have appeared in Ref. [6] (for k = 3) and recently in Ref. [8], with perfect agreement [19]. 

n 
F 

parafermions (but of course does not describe the total orbital angular momentum of the states). Those expressions 
now emerge from analyzing polynomials in z's. 



The corresponding result that generalizes 



for q 1 was given in Ref. [7] as the "truncated character" for 



It is possible to take the limit — > oo with m a fixed. 
In this limit, all n a — > oo, and by assuming that \q\ < 1, 



n 



i 



1 (1 - qm) ■ 



(19) 



Then the resulting generating function for the number of 
zero-energy wavefunctions with no upper limit on the de- 
gree in each variable is identical with the results of Refs. 
[14,15]. These references also explain how to shift the 
degree in q so that the subsequent limit N — > oo can be 
taken for each term of fixed degree in q. For the simplest 
case, one divides by q N l k ~ N ; and the limit counts the 
number of excitations at finite change in angular momen- 
tum of the edge of a very large drop of particles, similar 
to the results of Ref. [20] for the paired states (the k = 2 
case). In this limit the leading term in the generating 
function when N is divisible by k is q° and corresponds 
to the vacuum of the conformal field theory for the edge, 
as there are no quasiholes in the interior. More gener- 
ally, one may also take the and N — > oo limits with 
some quasiholes held fixed at the center of the drop, and 
arrange that the powers of q in general are the conformal 
weights of the states on the edge. For up to k quasiholes 
at the center, this was done in Refs. [14,15], and pro- 
duces the character of the chiral conformal field theory 
for each sector of edge states (there are k + 1 distinct 
sectors for the M = case treated here, but this number 
does depend on M [2]). 

To conclude, we now have a satisfactory understanding 
of the trial quasihole wavefunctions of the RR states. We 
expect that the approach can be generalized to obtain 
analogous results for the quasiholes of other quantum 
Hall phases that have trial states that are zero-energy 
eigenstates of short-range interactions, such as bilayer 
states, spin-polarized ground states of spin-1/2 particles, 
and the various spin-singlet clustered states for particles 
with spin 1/2 or 1 that have been introduced in the lit- 
erature; the results of Ref. [21] are relevant to the latter. 
We want to emphasize that ways of writing the quasihole 
wavefunctions other than the one presented here could 
still be of potential interest for the purpose of gaining 
useful physical insight, or for performing further calcula- 
tions of the properties of the quasiholes. 
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Note added: After submission of the original version 
of this paper, some additional facts have come to light 
which it may be useful to describe here for completeness. 
First, I have become aware that an earlier, difficult paper 
by Feigin and Loktev obtains some of the results herein, 
including the counting formula, eq. (18). In addition, 
there is a connection with other recent work on the zero- 
energy states for the k + 1-body Hamiltonian. Several 
authors [23-25] (some of whom concentrated on k = 2 
only) have discovered that for the system on a torus or 
cylinder, as the radius is taken to zero, an effective de- 
scription emerges that is combinatorially simple. [This 
builds on earlier corresponding results for the Laughlin 
Abelian (k = 1) states, and others [26].] For the cylin- 
der (which topologically is the same as the sphere), we 
will consider the single-particle orbitals (corresponding 
to single-particle L z eigenstates z m on the sphere) to 
be labeled by m + 1 = 1, . . . , + 1, arranged in se- 
quence on a line. Then the recipe is that, for M = 0, 
the Hamiltonian stipulates that in a zero-energy state, 
any two neighboring orbitals can be occupied by a total 
of not more than k particles, which can be distributed in 
any way between the two orbitals. (For the M = 1 case, 
the first one relevant to fermions, this becomes no more 
than k in any group of k + 2 neighboring orbitals, with 
of course not more than one in each orbital because of 
the Pauli exclusion principle.) For the torus, the same 
Hamiltonian applies for orbitals with periodic boundary 
conditions, + 1 = 1, with "neighboring orbitals" in- 
terpreted cyclically. Clearly, either of these specifies a 
combinatorial problem in which one would like to count 
the total number of such states for each N, N^, k, and 
one would hope that the result for the sphere (cylinder) 
agrees with the formula, eq. (18), obtained here. 

This type of problem is connected with a long history 
of results in combinatorics, for which see Andrews' book 
[27]. The problem for bosons in the cylinder/sphere case 
can be related to counting partitions, by identifying the 
length of each row of (the Ferrers- Young diagram of) a 
partition with the orbital number m + 1 = 1, 2, of 
a corresponding particle. As partitions are defined to 
have rows of (not strictly) decreasing lengths, each par- 
tition corresponds to a unique many-particle basis state 
for a system of bosons, the number N of bosons being 



G 



the number of rows of non-zero length in the partition. 
(Similarly, for fcrmions there is a correspondence with 
partitions with all rows having distinct lengths. There is 
a correspondence between the two problems, obtained by 
adding TV — i to the length of the ith row of a partition 
for the boson case to obtain the corresponding partition 
for fermions.) Notice that the total number of boxes Af 
(the number being partitioned) is related to the total L z 
by NN^/2 + N - L z = Af for M = 0. Then in the limit 
N<f, — > oo, in which there is no upper limit on lengths 
of rows, the simplest case is the k = 1 case for bosons. 
The zero-energy states then correspond to partitions in 
which two neighboring rows have lengths differing by at 
least two. The number of these is described by a gen- 
erating function that can be written in two ways; this 
is the celebrated Rogers-Ramanujan identity, which is 
thus connected combinatorially with the quasiholcs of the 
Laughlin v = 1/2 state in the plane. The generalization 
to k > 1 is given by the Gordon- Andrews identities, see 
Ref. [27]. One side of each of these identities agrees with 
the Nj> — > oo limit of our result, studied earlier in [14,15]. 
(There are further results also, in which some number I, 
< I < k of quasiholes are fixed at the origin [14,15]. 
In the chain of orbitals, this corresponds to a restriction 
that the occupation of the first orbital is < k — I.) 

To make full contact with our results, we require results 
for partitions with the additional restriction of a maxi- 
mum length for the rows. For the k = 1 case, this "fini- 
tization" of the Rogers-Ramanujan identities was found 
in Ref. [28]. For k > 1, the number of partitions obey- 
ing the given rules, and the appropriate finitization of 
the Gordon- Andrews identities, was derived in detail by 
Warnaar [29] in a manner very close to the above inter- 
pretation as particles in orbitals. [The identities are also 
given in Ref. [22], in a "grand canonical" form in which 
one multiplies the results for fixed N by x N (where x is 
another indeterminate) and then sums over TV with 
fixed.] Here boundary conditions can be imposed at the 
two ends of the chain, that the occupation of the first 
orbital is < k — l, and that of the last is < k — V. For the 
quantum Hall states on the sphere, these correspond to 
fixing I quasiholes at the north pole, V at the south pole. 
One side of each of these identities (specialized to the 
case I = I' = 0) agrees with the formula (18) (times q N ). 
These results thus confirm that the counting produced 
by these simple one-dimensional Hamiltonians [26,23-25] 
reproduces that obtained from the two-dimensional anal- 
ysis here. 

For the torus, the periodic boundary condition on the 
one-dimensional system can also be handled. The ends 
of the open chain can be joined and the restriction of to- 
tal occupation < k for the pair of orbitals 1, N^, can be 
imposed by taking Warnaar's formulas, setting l + l' = k, 
and summing over / to obtain the total number of zero- 
energy states. It is an open problem to derive this using 
the wavefunctions on the torus in a fashion similar to 



that of the present paper. 
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Appendix A: The proof that a symmetric polyno- 
mial that is non-zero when z\ = Z2 = ■ ■ ■ = Zk, and 
zero when in addition Zk+i = Z\, vanishes quadratically 
{zk+i — z\) 2 (or faster) as Zk+\ — ► z\, goes as follows. 
We will work on the sphere [13], using the homogeneous 
or spinor particle coordinates Uj, Vi, in terms of which 
Zi = Vi/ui for each i. We will use auxiliary coordinates 
like the wi in eq. (2) to form coherent states; these cor- 
respond to spinors ai, Pi with wi = Pi/ai (examples of 
such functions appear in the main text). In this form, 
any wavefunction is rotationally invariant if the rotation 
acts on all of the ai, Pi pairs as well as on all the Ui, 
Vi pairs. That is, the wavefunction is a sum of prod- 
ucts of factors of the forms u^j — v^j (or similarly with 
ai, Pi pairs in place of one or both particle spinors). 
This is because this singlet combination of two spinors 
is the only way to construct invariant functions out of 
the N spinors ttj, Vi and auxiliary spinors ai, Pi. For 
our purposes, the polynomials must be homogeneous of 
degree in each m, Vi (corresponding to each particle 
having angular momentum N$/2), and also symmetric 
under permutations of the particles. Then a wavefunc- 
tion can be analyzed into terms according to the angular 
momentum L of some set of k + 1 particles, schemati- 
cally * = J2l,l z ^L,L z {k + 1)^l,-lA n -k-l), where 
^L,L z (k + 1) is a function of the first k + 1 coordinates 
only, and ^l-l z {N — k—l) is a function of the re- 
maining particle and all the auxiliary coordinates. These 
are combined so that the state ^ is a singlet. If the state 
^L,L z (k+l) for L = Z/ max = (k+l)N,p/2 is nonzero, then 
in particular its component with L z = L is nonzero, and 
this is represented by Yii u f* > so that it is nonzero when 
all Vi = (i.e. z\ = ■ ■ ■ = Zk+i = 0). Hence we require 
these L = L max components to vanish. The rate at which 
the function vanishes as z/c+i z± = Z2 = ■ ■ ■ = Zk = 
is determined by the largest angular momentum value L 
for which the component "J l.l z (k+ 1) is nonzero. By an- 
alyzing the totally symmetric states of k+ 1 spins of mag- 
nitude A^/2, one finds that the largest values are L max , 
L max — 2, ... (the multiplicities are one for these first two 
terms, but larger than one for subsequent terms). Thus 
the first subleading term contains factors (UiVj — ViUj) 2 , 
or (zi — Zj) 2 , which is what we wanted to prove. This 
proof uses the sphere, but can presumably be repeated 
directly for the plane by reference to the center of mass 
as well as the total degree. However, this is unnecessary, 
as the counting of states for the plane subject to an up- 
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per limit N$ on the degree in each variable is always the 
same as for the sphere. 

From this result, it follows that if coordinates Zk+i, 
Zk+2, ■ ■ ■ , Z2k approach z\ = • • • = Zk, then the wave- 
function must vanish quadratically in each of the second 
k coordinates, and so if Zk+i = • • • = z 2 fc, it vanishes as 
(zi — Zk+i) 2k as Zk+i — ► zi, as was also used in the form 
of the residues. 

Appendix B: Here we describe some alternative 
quasiholc wavcfunctions, which are natural generaliza- 
tions of the ground state and quasiholc wavefunctions 
I 

*;v,j\^,fc(zi, ■■-,z N ;w 1 ,..., w n ) 



In the products Y\ (a> (for a = 1, . . . , k), i, j run over 

T,0:/3< a N + 1 ' Y,t3:t3<a N l3> and als0 1 runS 0ver 

E/3:/3< Q "/3 + 1 > ■ ■ ■ ■> E/3 : /3<a ™& ■ To ensure that each 
function is of the same degree in all coordinates Zi, we 
impose 2N a + h a = (2N + n) / k for all a (which are equiv- 
alent to the earlier conditions). Note that each function 
need not be of the same degree in all the w's, and one 
should not symmetrize over the w's. It is easy to see that 
when any k + 1 z's are equal, the function vanishes, be- 
cause in any term in the symmetrization over the z's, at 
least two must be of the same "type" a. In the limit as 
Z( P -i)k+i, Z( p -!) k +k Zp (forp = 1, . . . , N k = m fe ), 
the function vanishes if a further set of k coordinates co- 
incide with one another. Thus each of them reduces to 
the form of a non- vanishing residue with F — N — kNk, 
but in which $f,n' ,k' is of the form (20) with the param- 
eters as given, and which, unlike the residues above, still 
depends on the first n — fik w;'s as well as F of the z's. 
In this function ^fw v , the partition of N is replaced 
by a similar partition of F, with F a = N a — Nk (a = 1, 
. . . , k'), while the partition h a of n is unchanged, except 
that hk must be dropped. 

The fact that ^fn' y nere depends on some of the 

' 0' 

same w's as appear in front of it in the residue makes it 
nontrivial to verify that each nonvanishing residue can 
be obtained from some (linear combination of) functions 
of this form. One would like to expand the non- vanishing 
residues that result from eq. (20) in linearly-independent 
polynomials in w's, in order to find the dimension of 
the space spanned by the coherent states. The clustered 
part of the residue produces a linearly-independent set 
of (i N - F )/ k + n \ symmetric polynomials in the w's, while 

the remaining factor \I> f ,n' ,k' would be assumed by in- 
duction to be spanned by some linearly-independent set 
of #fn' y polynomials in the w's. However, we were 
not able to establish that the total number of indepen- 



in the form given by Cappelli et al. [3] (rather than 
those of Ref. [2]). In these partitions of N and n are 
again used, but are more conveniently labeled in a dif- 
ferent, but related way: N (resp., n) is partitioned into 
Ni > N 2 > ■ ■ ■ > N k (resp., hi < n 2 < ■ ■ ■ < h k ) with 
^' =1 JV a = N (resp., Yla=i ™ a = n )- ^ ne relation to 
the numbers used in the main text is m a = N a — N a+ i 
(Nk+i = 0), n a = *52p-p< a hf3. The wavefunctions are 
functions of the z's and only one set wi, . . . , w„ of quasi- 
hole coordinates. The coherent state wavefunctions arc 
given by: 



dent polynomials in w's in the products of these two sets 
is given by eq. (6). The reason is that the set of poly- 
nomials in w's obtained by multiplying one member of 
one of these two sets of polynomials by a member of the 
other may not be linearly independent, even though both 
the original sets were. (This may be formulated math- 
ematically as asking whether the set that spans the ex- 
pansion of ^ p n' y is linearly-independent over the alge- 
bra of symmetric polynomials, not just over the complex 
numbers.) We were not able to establish the linear inde- 
pendence of the resulting set in general, or that the num- 
ber of independent functions equals the required number, 
though in special cases of small systems it does [19]. If 
N' k' could be expanded in S„-harmonic polynomials 
in w's, then the argument would go through. As usual, 
S„ denotes the symmetric (or permutation) group on n 
objects. The S„-harmonic polynomials are a linearly- 
independent set of n! polynomials in the n variables w; 
which, when multiplied by arbitrary symmetric polyno- 
mials, span the full space of polynomials in w's (thus, 
they form a basis for the latter vector space over the al- 
gebra of symmetric polynomials) [30]. Though, in the 
cases we checked, there are enough S„-harmonic polyno- 
mials in each degree for this to work (and one must note 
that permutations of the w;'s in the wavefunctions (20) 
leave the functions of z's invariant, so only one member 
of each orbit under S„ in the S„-harmonic polynomials 
can be used), we do not know if the functions can al- 
ways be expanded in this way. Thus, while we suspect 
that the quasihole wavefunctions in eq. (20) do span the 
space VN.N4,,k for all N, N$, k, we have been unable to 
prove this so far. 
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